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Abstract 

Charmonia suppression is one of the main signatures for the formation of a deconfined medium. 
However, also nuclear effects, not related to the production of a hot medium, can affect the i/if/ 
yield. The determination, from the study of p-A collisions, of the i/ij/ behaviour in nuclear mat- 
ter is, therefore, extremely important to correctly quantify the amount of charmonia suppression 
induced by the deconfined medium. In this paper the new NA60 results collected at 158 GeV in- 
cident energy, i.e. under the same kinematical conditions as the In-In (NA60) and Pb-Pb (NA50) 
data, are presented and compared with p-A measurements from other fixed target experiments. 
Results on A-A collisions are also reviewed taking into account the new available information 
on the influence of cold nuclear matter on the J /iff production yield. Finally, results on the J/i/f 
polarization are shown for p-A and A-A collisions. 



1. Introduction 

J /iff suppression in A-A collisions is considered to be one of the main signatures for the phase 
transition from an hadronic to a deconfined medium lU. However, the J/i/f production can be 
affected also by cold nuclear matter effects, as, for example, the final state absorption of the cc 
pair in the nuclear medium, the parton shadowing or the initial and final state energy loss. The 
amplitude of these effects, not related to the formation of a deconfined medium, is determined 
from p-A collisions and then extrapolated to A-A interactions to be compared with the measured 
J /tff yield. 

At SPS energies, cold nuclear matter effects are usually parameterized, in the frame of the 
Glauber model, by fitting the A-dependence of the J /iff production cross section per nucleon- 
nucleon collision. These nuclear effects are quantified through the value of a*j? extracted, up to 

now, from the NA50 data at 400/450 GeV (crj£ = 4.2 + 0.5 mb) |Q]. Another way to param- 
eterize cold nuclear matter effects is based on a fit to the A-dependence of the J /iff production 
cross section using a A" power law. A value of a different from unity indicates how much the 
J/i/f yield is modified by the nuclear medium. It should be noted that both a and crj^ are effective 
quantities, since they represent the amount of cold nuclear matter effects reducing the J /iff yield, 
but they do not allow to disentangle the different contributions (e.g. shadowing, nuclear absorp- 
tion) playing a role in this reduction. Nuclear effects, evaluated from p-A data at 400/450 GeV, 
are then extrapolated to A-A collisions at a lower energy (158A GeV) assuming in both cases a 
scaling with L, the mean thickness of nuclear matter seen by the cc pair in its way out through the 
nucleus, and imposing crjf to be energy independent. The expected J /iff yield is then compared 
with the measured one, as a function of centrality. Following the aforementioned approach, a 
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further suppression (the so called "anomalous" J/if/ suppression), exceeding cold nuclear matter 
effects, is indeed observed at SPS in Pb-Pb [33 and In-In fl] collisions. It is clear that to correctly 
quantify the amount of suppression due to the formation of a hot medium, the J/iff behaviour in 
the normal nuclear matter must be determined with high precision. 

Several results on J/0- production in p-A collisions at fixed target experiments are available, 
covering different kinematic and energy domains. HERA-B at HERA, for example, has studied 
p-Cu, Ti, W reactions at 920 GeV m, E866 at FNAL p-Be, Fe, W collisions at 800 GeV H, 
NA50 at SPS proton induced collisions on several nuclei (Be, Al, Cu, Ag, W, Pb) at 400/450 
GeV Q] and finally NA3 at SPS p-H 2 , Pt collisions at 200 GeV 0]. It is important to note that 
none of these existing results has been obtained under the same conditions, in terms of energy and 
kinematic domain, as the SPS A-A collisions collected at 158A GeV. This implies that several 
assumptions (as the energy independence of cr^v ) had to be done, in order to extrapolate cold 
nuclear matter effects from p-A to A-A collisions. 

2. NA60 p-A results 

For the first time the NA60 experiment has studied the J/0 production in p-A collisions at 
158 GeV 081], in order to provide a reference collected under the same kinematic and energy 
conditions as the A-A data. 

The experimental apparatus of NA60 is based on a traditional muon spectrometer coupled 
with a vertex telescope made of Si pixel planes, close to the target region. The matching, based on 
the tracks coordinates and momentum between the tracks reconstructed in the two spectrometers, 
allows an accurate measurement of the muon kinematics and therefore an improvement in the 
quality of the results with respect to previous experiments. For details on the apparatus see 
The target system of NA60, during the p-A data taking, was based on 7 different targets (Be, Al, 
Cu, In, W, Pb, U) simultanously exposed to the beam. NA60 has also collected data at 400 GeV, 
with the same experimental setup as the one of the 158 GeV data taking period, to be used as a 
cross-check of the NA50 results taken at the same energy. 

2.1. J /0 kinematical distributions 

The investigation of the J/0 kinematical distributions may help to obtain further insights in 
the understanding of the charmonium production and initial state effects. In Fig [T| the rapidity 
(Jcm) distributions corresponding to p-A collisions at 158 GeV and 400 GeV are shown. A 1-D 
acceptance correction has been performed assuming realistic transverse momentum distributions 
and a flat cos6>cs shape (see Sec|4] for this definition). Within the narrow rapidity coverage 
of NA60 (one rapidity unit), the distributions can be reproduced by a Gaussian function. A 
simultaneous fit to the distributions corresponding to the different p-A collisions indicates that, 
at 158 GeV, the data can be described by a gaussian centered at midrapidity (yu y = 0.05 + 0.05) 
with a width <x y = 0.5 1 + 0.02. At 400 GeV, because of the wider distributions, the peak position 
is less constrained. In this case the fit is performed imposing the slightly negative mean value 
{Hy = -0.2) measured by NA50 with a high statistics data set at the same energy [23- As expected, 
because of the higher energy, a larger rapidity width (<x v = 0.8 1 + 0.03) is obtained, in agreement, 
within errors with the NA50 value (<x v ~ 0.85). Unfortunately, the precision of these results do 
not allow to investigate the negative rapidity shift observed by NA50, neither to confirm the 
smaller A-dependent shift of the center of the Xp (and therefore y) distributions observed by 
HERA-B H. 
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Figure 1: acceptance corrected rapidity distributions at 158 GeV (left) and 400 GeV (right). 



The J /if/ transverse momentum (pj) distributions have also been studied. The spectra are 
corrected with an 1-D acceptance, obtained assuming a flat pj distribution and realistic y and 
cosOcs shapes. A broadening of the pj distributions as a function of A is measured at both 
158 and 400 GeV. This confirms previous observations from other experiments f^, 0, OH EH 



>pp 



12 , HI, Hil] and is usually interpreted as initial state multiple scattering of the incoming gluons 
(Cronin effect) before the hard scattering which will produce the cc pair. The average < p\> 
dependence on the mass number can be parameterized with < p\ >=< pj > pp +p(A^ 3 - 1) [5] 
where A 1 ^ 3 - 1 is roughly proportional to the length L. As shown in Figf2](left) the < p\ 
values are compatible with a linear growth with the square of the center of mass energy s. On 
the other hand, the slope of the parameterization p shows an almost energy independent pattern, 
apart from a hint of decrease at low energy, corresponding to the NA60 158 GeV result. 

2.2. Cold nuclear matter effects 

New NA60 results on the nuclear effects affecting the J/if/ at 158 GeV have been presented 
at this conference [8]. Results, integrated over the full pj range are given in the rapidity region, 
covered by all the subtargets, corresponding to 0.28 < ycM < 0.78. Nuclear effects are evaluated 
comparing the cross section ratio, cr^ / 'crj^, between the target with mass number A and the 
lightest one (Be). By computing this ratio, the beam luminosity factors cancel out, apart from a 
small beam attenuation factor. However, since the sub-targets see the vertex telescope under a 
slightly different angle, the track reconstruction efficiencies do not completely cancel out. There- 
fore an accurate evaluation of such quantities and their time evolution has been performed target 
by target, with high granularity down to the single pixel level, and on a run-per-run basis. 

J /if/ cross-section ratios are shown in Figf2](right) as a function of L. In the same figure, also 
results from a similar analysis performed on the 400 GeV data sample are shown. In this case, the 
results refer to the kinematical region -0.17 < ycM < 0.33, corresponding to the same rapidity 
range, in the laboratory frame, as the 158 GeV one. Systematic errors include uncertainties on the 
target thickness, on the y distribution used in the acceptance calculation and on the reconstruction 
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Figure 2: Left: energy dependence of the < pl^, > pp (top) and p (bottom) parameters. Right: cross section ratios for 
p-A collisions at 158 GeV (circles) and 400 GeV(squares), as a function of L. 



efficiency. Only the fraction of systematic error, not in common to all the points, is shown, since 
it is the one affecting the evaluation of the nuclear effects. 

Performing a Glauber fit to the data, or using the a parameterization, the following results are 
obtained for the new NA60 p-A data: cr 3 J* = 7.6 ± 0.7 (stat.) + 0.6 (syst.) mb (a = 0.882 + 0.009 
± 0.008) at 158 GeV and cr^ s = 4.3 + 0.8 (stat.) + 0.6 (syst.) mb (a = 0.927 + 0.013 ± 0.009) at 
400 GeV. It has to be stressed that the crj^ result at 400 GeV is smaller with respect to the one 
extracted from the 158 GeV, pointing to an energy dependence of this quantity. Furthermore, the 
obtained value is in very good agreement with a previous result obtained by NA50 at the same 
energy 0]. 

The comparison with results from previous experiments (HERA-B, E866 and NA50) can be 
done in terms of the extracted a [8] or c?jf values, as a function of the Feyman x variable (xp). As 
shown in Figj3](left), where all the available measurements are compared, a strong dependence 
of <j^ s on xp and on the beam energy is clearly visible. Cold nuclear matter effects are stronger 
at high xp and for a fixed xp value their importance increases while lowering the beam energy. 
As shown in the figure, the new NA60 results at 400 GeV confirm the NA50 values obtained at 
a similar energy. On the other hand, the 158 GeV data seem to point to higher crj^ values, with 
a hint for an increase in the narrow explored xp region. It is also worthwhile to note that older 
NA3 results 17[] on J/ip production are in partial contradiction with these observations, showing 
lower crj^ values, rather similar to those obtained with the higher energy data samples. 

Although not shown here, the results do not present any scaling with X2, the fraction of 
nucleon momentum carried by the target colliding partons. Since X2 is the variable driving the 
shadowing effects, the break of this scaling confirms the fact that parton shadowing can not be 
the only mechanism invoked to describe the data. It is, therefore, clear that the interpretation 
of the kinematical dependence of the cold nuclear matter effects is extremely complicated since 
there are many competing mechanisms affecting the ]/tf/ production and propagation into the 
nuclei. A theoretical description of the cold nuclear matter effects over the full kinematic range 
is extremely complex and not yet available ll5Ulttll7ll . 



As discussed, the a or the crj^ values are effective quantities used to describe the convolution 
of all the cold nuclear effects not related to the formation of a deconfined medium. At this 
conference, a first attempt to take explicitly into account the parton shadowing contribution to the 
NA60 data has been proposed. It must be noted that at SPS energies, the charmonia production 
explores a range of x (the fraction of the nucleon momentum carried by the parton) corresponding 
to the antishadowing region, where parton densities in the nuclei are enhanced with respect to 
those of the free nucleons. The antishadowing, evaluated with the EKS98[18] parameterization 
for the nuclear modification of the parton distribution functions (PDF), causes an enhancement 
of the charmonium production cross section per nucleon-nucleon collision. Therefore, if this 
initial-state contribution is explicitly taken into account, a larger cr^ s value is needed to describe 
the measured data (cr^(158GeV) = 9.3 + 0.7 ± 0.7 mb and cr^(400GeV) = 6.0 + 0.9 + 0.7 
mb). Results depend on the adopted parameterization of the PDF nuclear modifications. Slightly 
higher crjf values (~ 5 - 10%), for example, are obtained if the EPS08 parameterization is 
used. 
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Figure 3: Left: Compilation of er . values versus xp. Right: J/i/i suppression pattern in In-In (circles) and Pb-Pb 
(triangles). Boxes around the points correspond to the correlated systematic errors, while the filled box on the right is 
the uncertainty on the absolute normalization of the In-In points. A further global 12% error, not shown, is due to the 
uncertainty on crj^ . 



3. Anomalous J /if/ suppression in In-In and Pb-Pb collisions 

As previously discussed, the cold nuclear matter reference used up to now for A-A collisions 
by NA50 [3] and NA60 Jl] was based on the value obtained from p-A collisions at 400/450 
GeV, assuming crj^ s to be independent of the incident energy and ycM range. As shown in Figj3] 
(left), this is not the case, since the cold nuclear matter effects strongly depend on the energy and 
kinematical domain. Therefore the evaluation of the size of the effects affecting A-A collisions at 
158A GeV must be based on the new p-A NA60 results at the same energy. The expected shape 
for the i/if/ distribution in cold nuclear matter is computed within a Glauber model, assuming 
the J /if/ production as a function of centrality (determined measuring the energy, Ezdc, released 
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in a Zero Degree Calorimeter) dNj^/dEzDC to scale with the number of binary collisions. The 

nuclear effects affecting the J/if/ are then implemented assuming a^J? = 7.6 + 0.7 ± 0.6 mb, as 
discussed in the previous section. This reference curve is then compared with the measured In-In 
dN™ e ° s /cIEzdc distribution, following the procedure described in fl]. As expected, since the 

value directly obtained at 158 GeV is higher than the one extracted from p-A data at 400/450 
GeV, the amount of anomalous suppression will be smaller with respect to the previous estimates. 
Furthermore, in the procedure followed up to now to build the reference curve, shadowing effects 
were not explicitly taken into account. While in p-A data only partons in the target are affected 
by shadowing, in A-A collisions also the projectile is involved and this further contribution has 
to be included when extrapolating the cold nuclear effects from p-A to A-A B2011 . It turns out 
that if the shadowing is neglected in this extrapolation, a small bias is introduced, resulting in a 
~5% artificial contribution to the suppression of the J/if/ yield (if the EKS98 parameterization is 
used for the shadowing evaluation). Therefore, if shadowing is properly taken into account in the 
extrapolation from p-A to A-A, a further reduction of the i/if/ suppression is observed. This is 
shown in Fig[3] (right), where the final ratio of the measured and expected ]/if/ yield is presented 
as a function of the number of participants (N parl ) for both In-In and Pb-Pb data. Even within 
this new definition of the reference curve, for very central Pb-Pb collisions (N part > 200) an 
anomalous ]/if/ suppression, of the order of 20-30%, is still visible. 



4. J /iff polarization 

The measurement of the J/if/ polarization is an important tool to clarify the quarkonium 
production mechanism, since different theoretical models predict different degrees of polar- 
ization 12 ill. For example, polarization is expected to be sensitive to the spin states of the cc 
pair [22], therefore its measurement can shed some light on the color-singlet and color-octet 
contributions to the production process. Contrary to other observables, as the differential cross- 
sections, in the polarization predictions uncertainties related to the theoretical inputs cancel out, 
therefore providing well-defined expected values. However, quarkonium polarization has always 
represented a puzzle, since, at collider experiments, the theoretical models, which were able 
to predict the observed J /if/ production cross section, clearly failed in the description of the pj 
dependence of the polarization II23I1 . 

Experimentally, the J /if/ polarization is measured from the full angular distribution of the 
quarkonium decay products: 

dN 2 v 2 

— — = 1 + Acos 6 + /j.sin26cos(f> + — sin'0cos2(f> (1) 

dcos8d<p 2 

where the polar angle is determined by the direction of the positive muon, in the quarkonium 
rest frame, and a chosen z polarization axis (defining the decay plane). <p is the azimuthal angle 
between the reaction plane, containing the colliding hadrons and the decay plane. The A parame- 
ter is traditionally called "polarization": zero value for A indicates that the 3/if/ are not polarized, 
while A-l(-l) means transverse (longitudinal) polarization. Non zero values of // and v indicate 
an azimuthal anisotropy of the distributions. A, fi and v depend on the chosen definition of the 
z axis. It is, therefore, useful, to determine these parameters in more than one reference frame. 
However, it can be noted that the knowledge of the full set of coefficients allows to analytically 
calculate their values in other frames, through appropriate transformations and knowing the J/if/ 
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kinematics |24]. Of course, this is not the case if only the A parameter is measured. In the litera- 
ture the most commonly used frames are the Collins-Soper (CS) one, where the z axis is parallel 
to the bisector of the angle between the projectile and target directions in the J/i/r rest frame, and 
the helicity (HE) one, where the z axis coincides with the quarkonium direction in the target- 
projectile center of mass frame. It is important to note that because of the A, fi and v dependence 
on a chosen reference frame, polarization results from different experiments can be compared 
only if the same frame is adopted. Most previous polarization analyses up to now were limited 
to the choice of only one specific reference frame and were restricted to the measurement of the 
polar angle distribution. On the contrary, recent results from HERA-B 12511 and NA60 [8J have 
been obtained measuring for the first time the full angular distribution of the decay products of 
charmonia and comparing different frames. 

NA60 results are obtained applying a 1-D acceptance correction, assuming realistic vcm and 
pj distributions and a flat cosO and <p spectra. The measured polarization values present a de- 
pendence on the J/ifr transverse momentum, as shown in Figj4] (left) for the helicity frame. In 
particular, HERA-B A results indicate an increase of the i/tfr polarization, moving from slightly 
negative values (corresponding to longitudinal polarization) at low pj to values close to zero 
(absence of polarization) at higher pj. The first polarization results obtained in p-A collisions at 
158 GeV and 400 GeV by the NA60 experiment confirm, within errors, the observed trend. 

In Figj4](left) results on the the azimuthal parameter v are also shown, always in the helicity 
frame. In this case no clear dependence on the i/ip kinematics is visible. Also the ft parameter 
is consistent with zero everywhere. Although not shown here, both HERA-B and NA60 results 
in the Collins-Soper reference frame give A values which tend to be more negative and larger in 
absolute value with respect to the ones measured in the helicity frame. This is an confirmation 
of the fact that the polarization parameters have a clear dependence on the chosen frame. 
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Figure 4: Left: ]/if/ polarization parameters A (top) and v (bottom) as a function of pr in the helicity frame for p-A 
collisions from HERA-B and NA60. Because of a slighty different definition of the v parameter, with respect to HERA- 
B, the plotted NA60 v values have been divided by a factor 2 with respect to those obtained from Eq. [T] Right: A (a,c) 
and y (b,d) results for In-In collisions versus pj and centrality respectively. The boxes around the points represent the 
total error affecting the measurement, assuming, for the NA60 points, a conservative preliminary 0.1 systematic error. 



The HERA-B and NA60 new polarization results clearly will shed some light on the polariza- 
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tion issue, in the attempt of describing all the available polarization measurement in a common 
scenario. 

NA60 has also provided, for the first time, results on the full angular distribution of the i/ift 
decay products in In-In collisions. The A values, in the helicity frame, do not show a dependence 
on the J /if/ pi as shown in Fig3] (a) and they are almost consistent with zero everywhere. Also 
the azimuthal coefficient v is rather small, pointing to a positive value only for pr ~ 1 GeV/c, as 
shown in FigS](b). The A and v values do not present any dependence even on the centrality of 
the collision, as shown in Fig|4](c,d). The pattern of these results is confirmed also if the Collins- 
Soper frame is adopted. In principle, it may be expected that the formation of a hot medium may 
affect the J /if/ polarization, as proposed in ll26ll . Therefore, quantitative predictions are needed, 
in order to clarify the observed results. 

5. Conclusions 

New NA60 results obtained in p-A collisions at 158 and 400 GeV have been presented and 
compared with the already existing measurements from other fixed target experiments. Cold 
nuclear matter effects, determined from p-A interactions, exhibit a rather strong energy and kine- 
matical dependence. Therefore the use of the new p-A results at 158A GeV, the same energy 
as the A-A data, allows a more precise definition of the expected cold nuclear matter effects in 
nuclear collisions. Apart from very central Pb-Pb collisions, where the anomalous suppression 
is still sizeable, a smaller effect with respect to previous estimates is observed. Results on J /if/ 
polar and azimuthal parameters have also been presented for both p-A and A-A collisions. 
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